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It is shown that the homogeneous coohng state (HCS) for a heavy impurity particle in a granular 
fiuid supports two distinct phases. The order parameter (j>a is the mean square velocity of the 
impurity particle relative to that of a fluid particle, and the control parameter ^* is the fluid cooling 
rate relative to the impurity collision rate. For ^* < 1 there is a "normal" phase for which (f)^ 
scales as the fluid/impurity mass ratio, just as for a system with elastic collisions. For 5* > 1 an 
"ordered" phase occurs in which (j)^ is flnite even for vanishingly small mass ratio, representing an 
extreme violation of energy equipartition. The phenomenon can be described in terms of a Landau- 
like free energy for a second order phase transition. The dynamics leading to the HCS is studied 
in detail using an asymptotic analysis of the Enskog-Lorentz kinetic equation near each phase and 
the critical domain. Critical slowing is observed with a divergent relaxation time at the critical 
point. The stationary velocity distributions are determined in each case, showing a crossover from 
Maxwellian in the normal phase to an exponential quartic function of the velocity that is sharply 
peaked about the non-zero (j)^ for the ordered phase. It is shown that the diffusion coefficient in the 
normal phase diverges at the critical point and remains so in the ordered phase. This is interpreted 
as a transition from diffusive to ballistic dynamics between the normal and ordered phases. 



I. INTRODUCTION 

A mixture of two mechanically different fluids rapidly approaches a common equilibrium state for times larger 
than a mean free time. This equilibrium state is characterized by a common temperature or, equivalently, mean 
square velocities for each type of particle that differ by their mass ratio according to the equipartition of energies. 
Recently, the corresponding state for a granular mixture was studied using a two component system of hard spheres 
with inelastic collisions [Q. Instead of the equilibrium state, the granular mixture attains a homogeneous cooling state 
(HCS) in which all time dependence occurs through a scaling of the particles velocities by their root mean square 
velocities. Although both components have a common cooling rate due to the inelastic collisions in the HCS, their 
granular temperatures are different. In terms of their mean square velocities, this implies a violation of the classical 
equipartition theorem. The extent of the violation depends on the mechanical differences of the particles (e.g., mass, 
diameter, coefficient of restitution), and is greatest when the differences are large. The quantitative predictions of the 
two temperatures from an Enskog-Lorentz kinetic theory have been confirmed by Monte Carlo simulations 1^. 

This effect also occurs for the simplest mixture of an impurity particle in a one component fluid. The impurity 
"equilibrates" to a common HCS with different temperatures for the impurity and fluid particles. The dynamics of 
an impurity particle of mass mo in a granular fluid with particles of mass m has been studied for the limiting case of 
m/mo ^ 1 The description was based on the Enskog-Lorentz kinetic equation for the impurity in a dense fluid 
and the fluid was taken to be in its homogeneous cooling state (HCS). As for the case of elastic collisions, the kinetic 
equation reduces to a simple Fokker-Planck equation in this limit with a velocity independent friction coefficient. 
The solution to this equation approaches an HCS for the impurity particle. As expected, the kinetic temperatures 
of the two types of particles (defined in terms of their mean square velocities) are always different although their 
cooling rates are the same. The ratio of impurity to fluid thermal velocities is not simply m/mo as for equipartition 
of energy, but has a more complex mass dependence, according to the mechanical properties of both particles and 
the degree of inelasticity in collisions. Nevertheless, the analysis requires that this mass dependence be such that the 
ratio of thermal velocities should vanish for m/mo ~* just as it would for equipartition. A single parameter the 
ratio of the cooling rate to the impurity-fluid particle collision rate, characterizes the domain for which the thermal 
velocity ratio vanishes, ^* < 1. The predictions of the Fokker-Planck equation in this domain (velocity distribution, 
temperature ratio, mean square displacement, diffusion coefficient) have been confirmed by both Monte Carlo and 
molecular dynamics simulation Q. As ^* ^ 1 the diffusion coefficient calculated from this Fokker-Planck equation 
diverges. 

The objective here is to put the analysis of reference |§1 in context by extending the discussion to ^* > 1. A 
preliminary report of this work has been given in reference H. It is found that there is a qualitative change in the 
state of the system at ^* = 1 that is analogous to a second order phase transition. The order parameter 0^ is the 
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ratio of thermal velocities with a conjugate field h proportional to the mass ratio. The parameter ^* is the analogue 
of the inverse temperature. The terminology "ordered" is used in analogy with magnetic systems where the ordered 
phase has a non-zero order parameter (magnetization) at zero external field. More precisely, the ordered phase here 
is associated with a broken symmetry or scaling Xh ^ Xcj)^ which applies for f * < 1 but does not hold for ^* > 1. For 
^* < 1 the fluid is "normal" with (p^ = at h — 0, a,s in the case of a system with elastic collisions. For ^* > 1 an 
"ordered" state with 0^ 7^ occurs at h — 0, representing an extreme breakdown of equipartition. Critical slowing 
and qualitative changes in the velocity distribution function for the impurity particle occur near the transition. The 
diffusion coefhcient diverges for ^* > 1 and can be understood as a transition from diffusive to ballistic motion. 

In the next section three characteristic frequencies are introduced: the cooling rate for the fluid particles, the cooling 
rate for the impurity, and the impurity-fluid collision rate. A simple estimate is obtained using a maximum entropy 
distribution to construct a phenomenological overvie w o f the HCS, its properties for m/mo <C 1 (or equivalently 
/i -C 1), and the phase transition analogy. In Section III the diffusion coefficient is calculated from its Grecn-Kubo 
representation using the leading term in a cumulant expansion of the velocity autocorrelation function The 
diffusion coefficient is expressed as a function of the order parameter (j)^ {^*,h), and for ^* < 1 the results of ||] are 
recovered. Otherwise, at the critical point and in the ordered phase, it is divergent. This divergence is interpreted 
by reconsideration of the Green-Kubo expression for finite times, showing a crossover from diffusive behavior in the 
normal phase to ballistic motion in the ordered phase. 

A more complete description is given in Section IV based on an exact asymptotic analysis of the Enskog-Lorentz 
kinetic equation for the impurity particle velocity distribution function. This distribution function is calculated in 
the critical domain showing a crossover from Maxwellian for ^* < 1 to an exponential quartic function of the velocity 
centered about a non-zero value for ^* > 1. The functional form of (psih,^,*) and associated critical properties 
are similar to those obtained in the phenomenological overview, with no qualitative differences. These results are 
summarized and discussed in the last section. 



II. PHENOMENOLOGICAL OVERVIEW 



Consider a fluid of hard, smooth, inelastic spheres of mass m, diameter a, and fluid-fluid particle coefficient of 
normal restitution a. In all of the following it is assumed that the fluid is in its HCS. Due to the inelastic collisions 
among particles the mean kinetic energy decreases as a function of time (referred to as "cooling" ) . An impurity 
particle of mass ttiq, diameter ctq, and impurity-fluid particle coefficient of restitution ao is inserted in the fluid at 
some initial time. There is energy transfer between the impurity and fluid particles due to collisions and subsequently 
a common HCS for the fluid and impurity is attained where all particles have the same cooling rate. In this section a 
phenomenological but accurate description of this process and the HCS is given to present the basic ideas in a simple 
physical context. 



A. Nonlinear friction coefficient 



The primary property of interest is the ratio of the mean square velocities for the impurity and fluid particles 

where the brackets denote an average over the initial state of the fluid plus impurity particle. This function measures 
the accommodation of the impurity particle to the fluid and will be referred to in the following as the order parameter. 
The cooling rates associated with the mean square velocities are defined by 

m^-dtHv'it)), Ut)^-dtHvl{t)). (2) 

For dimcnsionlcss units it is useful to define an average impurity-fluid particle collision rate 

, 8 _2 , , , , , 1 -l- ctQ m 

Mt) = o Vttct 50 w i ; = — ■ , 3 

6 2 m + mo 

where ct = (ct -t- cto) /2 is the average diameter, p is the fluid density, go is the pair correlation function for the impurity 
particle and a fluid particle at contact, and {v{t)) is the average speed of a fluid particle in the HCS. The parameter 
h has been introduced as a measure of the mass ratio. As a function of h this form for the collision frequency is the 
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same as that for elastic collisions characterizing the equilibration rate. A dimensionless equation for (j){t) now can be 
written in the form 



ds<l>={C-Co)cl^, (4) 
where the dimensionless cooling rates and dimensionless time have been introduced as 

r = -, ^0 = -: ds^Mt)dt. (5) 

To proceed it is necessary to calculate ^* and as functions of cj). As shown in Appendix these are related to 
averages over the pair distribution function for two fluid particles and for a fluid and the impurity particle, respectively. 
This is a formal result since the distribution functions are not known. As a phenomenological estimate therefore, 
these averages are performed using a maximum entropy ensemble parameterized by the true mean square velocities. 
The qualitative accuracy of this approximation is confirmed in Section IV. The results of Appendix |A] are 



A-s/2h go No- 
where g is the pair correlation function for two fluid particles at contact. This form for the cooling rate of the 
impurity is the same as that for elastic collisions and represents the equilibration rate. The new features of inelasticity 
are primarily described by f*, which is independent of (jj. The equation for (f){s) with these approximate forms for the 
cooling rates is 

[9,+r(0)-r]0 = w), (7) 

which results from the decomposition ^q(0) = 7*(0) — hn* {(j>) / (j> with the definitions 

r&)^{l+^f\ n(0) = (1 + 0)'/'. (8) 

This has the same form as would be obtained from a simple Langevin or Brownian motion model where 7*(0) is the 
"friction constant" or nonlinear impurity-fluid collision frequency and hn{(f)) is the noise amplitude. The solution to 
this equation is a function of time and the two parameters ^* and h. The stationary solutions (p^ (^*, h) are determined 
from 

4 = h . (9) 

This form shows most clearly the effect of competition between "friction" on the impurity particle and fluid cooling 
since 7* (0s) > £,* is required for positive, finite solutions. This generalizes the result obtained in reference [|| which 
is limited to ^* < 1 and ft, — > 0. It is easily verified that a unique positive solution to (||) exists for all positive 
^* and h and that it is linearly stable. The latter confirms that the HCS characterized by is approached for 
long times for a wide class of initial conditions. The time scale for formation of the HCS is discussed below. For 
elastic collisions {a = ao = 1) the solution is (j)^ = h/ {1 — h) = m/mo as required by equipartition. If only the 
impurity-fluid particle collisions are inelastic (i.e., a = 1, ^* = 0) a recent result of Martin and Piasecki is recovered 
§, 0s = V (1 - /i) = + ao)/ [2mo + m (1 - ao)]. 

More generally, Eq. (S) can be transformed into a cubic equation for 0^ whose physical solution gives 0s (^* , h) for 
arbitrary ^* and h. Figure |l| shows 0s {£,*, h) as a function of h for ^* = 0.9, 1.0, and 1.1. An instructive alternative 
form for the determination of 0s is 

C^CM- (10) 

The graphical solution to Eq. ([l^) is obtained in a plane y vs by finding the value of at which the constant 
y = £,* intercepts the curve y — as illustrated in Fig. ||. There is seen to be a qualitative difference between the 

solutions for ^* < 1 and ^* > 1 in the limit of small h. Since 7*(0) > 7*(0) = 1 and n(0) = 1 the asymptotic solution 
for r < 1 is 
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which agrees with [gj. The mechanism responsible for solutions with ^* > 1 is now clear. As 5* exceeds 7*(0) the 
nonlinear dependence of the friction coefficient on is activated to maintain positivity of 7*((^s) — Since 7*(<^) 
is a monotonically increasing function of (j), positivity is possible for any choice of ^* . In general this requires that </> 
must be finite even for ft, — > 0. This is possible if 7*(0s) — is of order h for small h or 4>^ — constant + 0{h). This 
nonlinear dependence of the friction coefhcient on (/)g provides the mechanism whereby the coupling of the impurity 
particle to the fluid can be enhanced for large cooling rates: the impurity-fluid collision frequency is increased by an 
increased mean square velocity of the impurity relative to that of the fluid. This is illustrated in Figs. |l| and ^ showing 
the qualitative difference between ^* < 1 and ^* > 1. The former admits ^ for h ^ whereas the latter requires 
^constant. In more detail, the asymptotic solution to (|^) is 

r (i-D^'/i, r < 1, 

cj^siCh)^! V2h, r = l, (12) 

The common domain of h —f and ^* ^ 1 can be obtained from ^ by the scaling h —>■ e^h, — > e^^, and 
1— ^* e (1—^*), for e <C 1. The result is the quadratic form 

ft « (1 - r)0s + (13) 

which has the solution 

« r - 1 + \/(e-i)' + 2ft. (14) 

At ft = this gives « ^* — 1 + |^* — 1|, illustrating again the qualitative difference between ^* < 1 and > 1. 
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FIG. 2. Plot of 5S(<^), Eq. (g), for h 
the curves with the horizontal lines ^* = 



: 10 (dotted line), h = 10 ^ (dashed line), and h — (solid line). The intercepts of 
0.9, 1, and 1.1 give the corresponding values of (j}^{£^*,h) (circles). 



Figures and ^ show that the fluid cooling rate relative to the impurity-fluid collision rate £,* is a control parameter 
distinguishing different dependencies of (j)^ on h for small h. This will be exploited in the next subsection, where f * = 1 
identifles a critical point. Since ^* cx (1 — a'^)/h, Eq. (H), the plots of (p^ at constant ^* require the change of both 
the fluid coefficient of restitution a and the mass ratio parameter h. It is instructive, however, to examine the mean 
square velocity ratio (j)^ and the mean energy ratio ea/e — {mo/m)4>g as functions of h at fixed a < 1. In that case, 
£* ~ diverges in the limit h ^ and so do cj)^ w + — 1 '--^ and eo/e ^ h~^. This is illustrated in Fig. 
y, where <j)^ and eo/e arc plotted versus h for a — 1, 0.99, and 0.95 (taking, for simplicity, ao ~ 1, a = ctq, g = go = !)■ 
The dotted lines for a — I represent equipartition for which (f)^ — )■ h/{l — h) and eo/e 1. For a < 1 there is a sharp 
deviation at sufficiently small h representing the crossover to the domain for which ^* > 1. 
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FIG. 3. Ratio of mean square velocities, 4>^, and of mean kinetic energies, eo/e, as functions of the mass ratio parameter h 



B. Representation as a phase transition 

Figure]^ and Eqs. ( |l^ ) and are reminiscent of the thermodynamics for magnetization as a function of an external 
magnetic field. Below some critical temperature the magnetization is finite at zero field, while above that temperature 
it vanishes at zero field. To pursue this analogy, consider (t>si^*T^) the order parameter (magnetization), h the 
conjugate field (magnetic field), and £,* as the control parameter (inverse temperature). The "equation of state" for 
the system is obtained from (ph by solving for h (c/j^,^*), 



,1/2 



,3/2 



(15) 
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A Helmholtz free energy can be defined in the usual way 

F r) = dx h (x, n - 4 - In (1 + - 2e* 



^1 + 4) 



1/2 



(16) 



Next, the Gibbs free energy is obtained from the Legendre transformation 



2 + 0s(r» 



1/2 



(17) 



The first and second derivatives of $ (^*, h) provide the order parameter (j>^, "entropy" S, "susceptibility" x, "expan- 
sion coefficient" ah, and "heat capacity" Ch- The results are 



2 + 4 



- 2 



(18) 



x(r,M = — iu^= = — _ \/2 (19) 



a/iC,/i) = ^ = X _ ^ (20) 

9e (i + <^s) 



a2$(e,/i) 



= H7I5^=X"«t (21) 



The values of these thermodynamic properties in the limit h = for ^* 7^ 1 follow directly from the asymptotic forms 
(H) for 4: 

$ (r, = 0) = I _ ^3 _ _ 2 In^*, I* t 1. ^^^^ 

f 0, r < 1, 



1, r < 1, 



xiC,h = 0)^j^::^{ (24) 



a,(r,/» = o) = ^ 2r, f*ti, (2^) 

C,.(r> = 0) = | (26) 

With the exception of thermodynamic variables vanish for ^* < 1 and are finite for ^* > 1. All are continuous 

at ^* = 1, except for ah which has a finite discontinuity. The susceptibility diverges as |C* ^ 1| ^ 0. Thus cither the 
discontinuity of ah or the divergence of x characterizes a second order phase transition at ^* = 1. Since the order 
parameter (p^ behaves qualitatively as that for a system with elastic collisions when ^* < 1, this will be referred to as 
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the "normal" phase. In contrast, since 4>s for ^* > 1 this will be called the "ordered" phase. The entropy function 
E(^*, /i) is plotted versus ^* for h = 10"^, 10~^, and in Fig. ||. The negative value of E at = and ^* > 1 is a 
measure of the degree of "order" in the ordered phase. The response functions x, ah, and Ch a-t h — 0, 10^^, and 
10^'^ are shown as fimctions of i^* in Fig. ||. 




V, 



FIG. 4. Entropy as a function of ^* for h — 10 ^ (dotted line), ft = 10 ^ (dashed line), and h = (solid line). 
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FIG. 5. Inverse susceptibility (x expansion coefficient {an), and heat capacity {Ch) as functions of ^* for h = 10 
(dotted line), h = 10~^ (dashed line), and h = (solid line). 
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Near the critical region {h <^1, |^* — 1| ^ 1), the free energy adopts the Landau- hke form 

HC,h) « i(l - at + - h4,, (27) 

which yields the critical equation of state (|l3|), as expected. It is easily verified that the free energy and the equation 
of state in the critical region satisfy the scaling relations 

$(A(r-l),A'^/i) = A''<i>(r-l,/i), UHC ~l),X''h) = X'--Ue -l,h) (28) 

with a = 2 and 6 = 3. These scaling relations suffice to determine the critical exponents Q S = a/(b — a) = 2, 
(3 = h — a = 1, and 7 = 2a — 6 = 1, while the critical exponent d==2 — 6 = — lis negative, indicating that Ch is 
continuous at the critical point. 

C. Critical dynamics 

If the ratio between the initial mean square velocities of the fluid and impurity particles is not that given by the 
solution to (||), there is an evolution to the HCS described by (|^) which can be written in the Ginzburg-Landau form 

ds(j) = -n{(t)) 4= (29) 

Here, $ (^* ,h;(j)^ is a variational free energy given by Eq. ( [TtI) with the order parameter (j) considered as an independent 
variable, and the kinetic coefficient is n{<j>)- The stationary solution occurs for a$ {Ch; (f>) /d(f> = 0, which is just Eq. 
(|). It follows directly from (|l^) and that $ {Ch;^) has the properties 

This shows that $ (^*, h; 0) is a Lyapunov function for the dynamics: it is bounded from below by the HCS solution 

and monotonically approaches this bound. Consequently, the HCS solution results in both phases for a wide class of 

homogeneous initial conditions and is stable. 

The free energy $ (^*, ft, = 0; 0) is shown in Fig. ^ for ^* = 0.9, 1, and 1.1. As expected, the minimum is located at 
= for ^* < 1 and at ^ for ^* > 1. For states near the HCS the evolution equation (|2^ ) can be linearized and 

a characteristic response time t* identified according to 

-dsln\$-^J=T*-' = {nx-\. (31) 

In the elastic limit r* is just the equilibration time (in terms of the number of impurity-fluid particle collisions) for the 
impurity particle to attain a mean kinetic energy equal to that of the fluid particles. Similarly, for inelastic collisions 
it is the time for the impurity particle to reach a cooling rate equal to that of the fluid. This characteristic time is 
a smooth function of h and £,* except in the limit ft. ^ where r* diverges at ^* — 1. This critical slowing follows 
directly from the fact that t* oc x- Otherwise, the relaxation times away from ^* = 1 are flnite and comparable for 
the normal and ordered states. Figure |^ shows the dependence of r*~^on ft for ^* = 0.9, 1, and 1.1. 



a$ (r,ft; 



< 0. 



(30) 
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III. DIFFUSION 



Diffusion of an impurity particle in the HCS has been described in general elsewhere P,p|,pl| . In this section the 
consequences for h in the two phases are explored. A generalized diffusion equation can be obtained by extending 
the familiar methods of linear response to the granular fluid, which for long wavelengths takes the form 

dsn*{r*,s) - D*{s)V^n*{r*,s) = 0. (32) 

Here n*(r*,s) is the dimensionless probability density to find the impurity particle at position r = r*£, where 

£ = {vq)V'^ /vc is an effective mean free path, and s is the dimensionless time of (^). The time dependent diffusion 
function -D*(s) is given exactly by a Green-Kubo expression 

D*{s) ^ -4^ ds'i^l {/) ■ v*)„ (33) 

where Vq — vo/-\/ (w^)s and the brackets denote an average over the dimensionless HCS ensemble. A phenomenological 
but accurate evaluation of the velocity autocorrelation function is given by its exact short time behavior 

{v^is')-v*,l^{v-^%e->\ (34) 



where L* is the dimensionless Liouville operator [cf. Appendix |A| . The dimensionless frequency is calculated in 
Appendix ^ using the same approximation as that for the cooling rates in Section |ll|, with the result 

^*D = l[r&s)~eo]- (36) 

For a fluid with elastic collisions the approximation (|3j) coincides with that obtained from the Enskog-Lorentz 
equation in the first Sonine approximation, and is known to be accurate even for moderately dense systems. It is 
assumed that a similar level of accuracy extends to the inelastic case as well The diffusion function D*{s) 

becomes 

..*(.) = ^(l-e--^). (37) 



The analysis of Appendix ^ shows that uj'^ > for all finite h. Thus for s ^ uj 



D 



D*{s)^D*^-L. (38) 
3w^ 

and ( ^^ becomes the usual diffusion equation with diffusion constant D* . The initial transient period is the expected 
"ageing" required for applicability of hydrodynamics (diffusion) . 
Consider now the behavior as h —^ 0. Using (^ and (^) cjj, can be expressed entirely in terms of and h 

/, -r x3/2 

= 2'^^ ■ (39) 



Using ( |12D this frequency behaves for h ^ as 

i(i-e), r < 1, 
^L(r,M- <! 5^1' r = i, (40) 
^,^,h, r > 1. 

In general, (^*,0) is finite below the critical point, but vanishes at and above the critical point for h = 0. Thus 
diffusion in the sense of (38) occurs at /i = only for ^* < 1. To understand the phenomenon for ^* > 1 note that for 
wjj = Eq. (M) becomes 
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= (41) 

To interpret this, take the second moment of ( ^2|) with respect to to relate D*[s) to the mean square displacement 
of the impurity particle, 

D*is)^ldA\r*is)-r*i0f)^. (42) 
Thus the mean square displacement behaves as 



(|r*(.)-r*(0)|^).-|^P' (43) 
This shows that the impurity is not diffusing but rather undergoing ballistic motion at its root mean square speed if 

r > 1- 

IV. ASYMPTOTIC KINETIC THEORY 



The analysis of Sections and [II is based on the plausible but uncontrolled estimate of the cooling rates for the fluid 
and impurity particles using a maximum entropy ensemble (Appendix]^. Furthermore, it is limited to a discussion 
of the order parameter and diffusion but does not address other properties such as the velocity distribution itself. In 
this section the results of Section || are recovered systematically and with additional detail from the Enskog-Lorentz 
kinetic equation for the impurity particle velocity distribution The features of interest here occur for h ^ 

so only an asymptotic representation of the kinetic theory is required. The fluid particle distribution is independent 
of h and its detailed form is not required for the analysis here. The asymptotic form of the Enskog-Lorentz equation 
for the impurity particle distribution /o(vo,t), as a functional of the fluid particle distribution, is the focus of this 
section. 

An expansion of the impurity-fluid particle collision operator in powers of h is straightforward, leading to the 
Kramers-Moyal representation |11[| . The leading terms of this expansion have been given in Appendix A of reference 

d 1 d'^ 

9t/o(vo,i) = [/ivo7(wo)/o(vo,i)] + -t; — ^ — {h'^ ["i(wo)^jj 

dvo 2 dvoidvoj 



+n2{vo) (^voiVoj - ^Sijvl^ 



Mvo,t)}+O{h'). (44) 



The friction ^(vq) and the noise functions ni{vQ), ri2(wo) are explicit averages over the fluid particle distribution given 
in Appendix M. The states of interest are functions only of the magnitude of vq. Consequently, it is possible to 
introduce a variable 

(45) 



" (vHt)) 

whose average value is the order parameter <j){s), where s is the dimensionless time variable defined in Eq. The 
distribution function for this variable is -P(0, s), defined by 

P{^,s) EE AnfoMvl^ = 27T{vHt)f'^'^'foivo,t). (46) 
a(p 

Then the Kramers-Moyal expansion becomes for P{(j>, s) 

a r ....... A 2 9 A , 4 9 



+0{h^). (47) 
Here ^* is the dimensionless cooling rate for the fluid introduced in (||). The functions 7* ((/>) , n* {(j))^ and 71-2(0) ^'"^ 



the dimensionless forms of 7(^0), "i(i'o), and n2(wo), respectively, given by Eqs. (A27) and (|B3) of the appendices 



L2[. They are functional of the distribution function f{v,t) of the fluid in the HCS and are normalized to have 
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7*(0) = 1, nt(0) = 3(?;3)/4(w2)('i;), and 71^(0) = 1. In addition, the derivative ■j*'{(f>) = dj* {(f)) / d(j) at (/) = is 
7*'(0) = {v~^){v'^) /5{v) . These functions can be ac cura tely estimated by assuming a maximum entropy ensemble for 
the fluid / and the results are given by Eqs. (B8)-( B10 ). According to these estimates, nl{0) ~ 1 and 7*'(0) ~ 3/10. 
Figure ^ shows the friction coefficient 'y*{(f>) and the noise coefficients ((/>), n2{(f>) according to this maximum entropy 
appro xim at ion f or /. Hereafter, all the plots of quantities defined in terms of those coefiicients will be made using 
Eqs. (^)- (IbT^) . 
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FIG. 8. Plot of the friction coefficient 7*(0) (solid line) and the noise functions n^{<j)) (dashed line) and n2{(j}) (dotted line). 

The Kramers-Moyal expansion is not well-ordered since the small parameter h also multiplies the highest (j) deriva- 
tive. A proper asymptotic result requires a scaling such that all higher terms in the series are exactly zero in the 
appropriate limit. The simplest case is the deterministic limit for which h = in (p7|). 

A. Deterministic limit 

If the formal limit h — is taken in (p7|), the equation becomes 

dsPoicI), s) = -^cj, [-C + l*m Poicl^, s), (48) 

where the subindex is used to denote quantities at h — 0. The solution to this equation for sharp initial conditions 

Po(0,s-O) -J(0-0o) is 

Po('^,s) = 5(0-0o(s)) (49) 

with 

[9,-t-7*(0o)-r]^o(s)=O (50) 

and 0q(s = 0) = 00- Thus the initial sharp distribution remains sharp and only its central value changes in time The 
latter defines the macroscopic dynamics for the average value <j>Q{s). As expected, it has the form (^ with a vanishing 
noise. The solution for more general initial conditions can be obtained as a superposition of the specific solution (|4S|). 
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The stationary solutions are obtained from (50) as the solution to 



(51) 



and 7* (00s) = C*- It is shown in Appendix ^ that 7*(0) — 1. Therefore, the solution 
< 1. In the case ^* > 1 the unique stable solution is determined from Yi't'os) = C with a 
non-zero value of 0os- Such solutions exist because "f*{4>) is a monotonically increasing function, i.e., j*{4>) > 7*(0), 
7*'(0) ^ Oi as proved in Appendix]^. These are the two phases discussed in Section now identified precisely from 
the Enskog-Lorentz kinetic equation. The details of the "equation of state" are different for this controlled analysis, 
but the qualitative features of states with = and </>s for h = are recovered exactly, as illustrated in Fig. H. 



The possibilities are (/jq^ - 
(pQg = is stable only if ^ 



Os 




FIG. 9. Plot of the order parameter in the deterministic limit, ^q^,, as a function of The dashed line is the maximum 
entropy estimate (^j, = — 1 of Sec. 



B. Effects of fluctuations 



A more complete description including fluctuations is obtained by a transformation of the form (j) = 4>q{s) + h^rj, 
where (jjois) is the average value of at /i = and 77 represents the fluctuations about this value. The power law 
of the scaling for the fluctuations is determined by the requirement that the distribution of fluctuations V {rj, s, h) = 
h^P{(l), s,h) ^ V (?7, s) which is independent of h. Inverting the result in terms of (j> gives the well-defined asymptotic 
behavior for small h ||ll|. Here, attention is limited to states near the stationary state t^Qg so the chosen scaling is 
(/) = -l- hPrj. The distribution function is no longer sharp, as in (^9|), but instead has a width proportional to h^. 
The choice of p is governed by the requirement that the Kramers-Moyal equation for the distribution of rj should 
truncate exactly for h — 0. The details are given in Appendix where the stationary solution in the normal phase 
is found to be 



3/30 



1/2 



-34>/2hJJ^ 



C < 1- 



(52) 
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and the width of the distribution is characterized by 



1-e 4(w2)(z;)i-e' 

Since (/)gg = in this phase, the order parameter is cj)^ — hrj^. This agrees with the resuh of reference where the 
distribution is recognized as a MaxweU-Boltzmann distribution for the velocity of the impurity particle, but with a 
different temperature than that of the fluid. In the present notation the impurity temperature identified from this 
Maxwellian is 

To = T^A^%. (54) 

where T is the granular temperature of the fluid. The phase transition is seen to occur with a diverging kinetic 
temperature for the impurity particle. If a maximum entropy distribution is assumed for the fluid, then the right side 
of ( p3| ) can be evaluated to get 0g = h/{\ — ^*), which agrees with the phenomenological theory of Section ||, Eq. 

®- 

In the ordered phase a qualitatively different distribution is obtained, as expected. It is now Gaussian in (quartic 
in velocity) and centered about a non-zero value 

P^(0) ^ ^=^_e-(*-^=)'/2^(^=)'\ r>l, (55) 



with 

B{<f) 



7*'(0) 



2 4 



(56) 



The function B {4>) is plotted in Fig. ^ The width of the distribution is A(/)s — [i?(0g)/i] , so that as /i — > the 
distribution becomes sharply peaked about the stationary value (j)^ — (jj^^, where is the stationary order parameter 
in the deterministic limit. At a fixed small value of h the (absolute) width A(f>s increases, but the relative width 
A(j)s/<j)s decreases, as ^* — 1 (and, consequently, 0g) increases. 
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0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 10. Plot of the function 5(0) defined in Eq. 



C. Critical domain 



The above results distinguish the cases of /i — > for ^* < 1 and for ^* > 1. A uniform description of the critical 
domain for small h and ^* w 1 can be obtained by noting that (p^ vanishes at the critical point from both phases, and 
scaling the Kramers-Moyal equation according to ^* — 1 = and (j) = h^^^r]. In addition a new time variable is 

defined by r = h^^^s. Then at ft- = the equation is 



drV{ll,T) 



d_ 



1 



-(5r; + 7*'(0V-- 3-2— 7, <(0) 



d 



drj 



The stationary distribution function is found to be 

Vs [t]) = Ct]'!'^ exp 



2B(0) 



7*'(0) 



(57) 



(58) 



where 5(0) = 2nt(0)/37*'(0) ~ 2.22 is the value at = of the function defined in Eq. (|6|). The scaled order 
parameter in this critical domain is then obtained from 



""^^^^ ^''^°^/„°^duuV2exp[-(.-2z)V2]' 



where z = S/2-f*'{0)y/ B{0). Its exphcit expression is 



(59) 



(60) 



where 
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r (3 + Az^)^2{z) + [Ky^z^) + X„i/4(z2) - K.^iiz^) ~ K_,^i{z^)] , z < 0, 

+ 6z^)^,{z) + 2z^ [Ir/,{z^) + I-7/i{z^) + h/d^^) + , z > 0, ^^'^ 

z <0, 



(62) 



^^)-i^3/.-^-^^ 
3/4[Z^) + 1-3/4 

I^{z) and Ku{z) being the modified Bessel functions of first and second kind, respectively. Tlie width of the distribution 



^ (z\-\ ^l/4(^') - ^3/4(^^), 



A77 = \^ rf — rj in the steady state can be obtained by taking moments in Eq. ( pT] ) as 



Ar;,(5) 



The asymptotic behaviors of ri^{5) and A?7s(5) are 



7*'(0) 



nm\5\-\ 5 
W(0), 



0, 
00, 



/273^(0)|(5|~\ (5^ -00, 
A77,(,5) ^ <( VA^, (5^0, 
/B(0), ^->«3, 



(63) 



(64) 



(65) 



where Ai = 2B(0) [r(5/4)/r(3/4)]^ ~ 2.43 and A2 = 3S(0)/2 - Ai ~ 0.90. In the limit 5 -> -00 we have Ayys/^, ™> 
•\/2/3, which is consistent with a Maxwell-Boltzmann distribution. In contrast, Ar/s/rJ^ — > when (5 — s- cx), so that the 
distribution is sharp around the order parameter in that limit. The dependence of the scaled order parameter rj^^S) 
on the scaled control parameter 6 is shown in Fig. where also the width of the distribution is plotted. 
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FIG. 11. Plot of rj^{5) (solid line) and Ai]s (dotted line). The dashed line represents the phenomenological approximation 
r^^{6) ^ 5 + VS^, Eq. (0). 

It is worth noting that the scaling relations (j) = h^^^rj, f * — 1 = h^^^S are successfully captured by the phenomeno- 
logical theory of Sec. |l|, except that there the scaling function is approximated by r]g{5) — > (5 + \/P~\~2, Eq. (pl|). 
While this function is not quantitatively correct, especially for d > (cf. Fig. ^l]), it is qualitatively consistent with 
the hmits in (|6|), the numerical coefficients being replaced by Ai ~ 2.43 -> 2, 7*'(0) ~ 0.3 7*'(0) = 1/2. 

Let us now go back to the unsealed variable (j). The corresponding distribution in the critical domain is 



Ps(0) oc (/)i/2exp 



1 



2B{Q)h 



C - 1 



and the equation of state is 



C - 1 



(66) 



(67) 



This result encompasses the normal and ordered phases, as well as the critical point. The normal phase in the critical 
domain is defined by /i <C 1 — ^* (i.e., S — oo), the ordered phase is recovered in the case /i ^ ^* — 1 (i.e., 6 — > oo), 
while the critical point corresponds to ^* = 1 ((5 = 0). Thus, the asymptotic behaviors (|6^) translate into 



hni{o)/{i~e), r <i, 
v^, r = 1, 
(e - i)/7*'(o), e>i. 



(68) 



In the normal phase, 



^ 1 — , so that the distribution ( pq ) becomes 

30(1 -D" 



Ps{cl)) OC 0^/2 exp 



2hnl{0) 



(69) 



which agrees with (p2[). In the ordered phase, however, the width of the distribution is much smaller than the average 
value = (^* — l)/7*'(0) [which is the solution to ^* = l*{4's) t;he critical region], so that the prefactor (j)^^'^ in 



(pq) can be replaced by a 



t1/2 



with the result 



(70) 



As expected, ( |70| ) agrees with 
is 



particularized to the critical region. Finally, at the critical point the distribution 



PM oc 0l/2e-0V2B(O)h^ 



(71) 



As anticipated from the behavior of Ar]s/r]s, there exists a crossover in the critical domain from the Maxwell- 
Boltzmann distribution (|69| ) to the sharp distribution (^oh through (71). Of course, the distribution function (|6^ ) is 
more general than the three limiting cases described by (|69|)-([7l|). To focus on the shape of the distribution function 
around its average value, define the normalized distribution P*{x) = 4>gPs{(f> = x^g) — rJ^Vsif] — ^Vs)- From Eq. (|5^ ) 
we have 



P*ix) oc a;^/^ exp 



- 2z 



^ <5/27*'(0)/B(0). 



(72) 



By construction, this distribution is normalized to (x) = 1, regardless of the value of the scaled control parameter S. 
The asymptotic forms of P* {x) are 



1 3x ~3x/2 
o [r(5/4)]^/^ ^l/2, 



[r(3/4)] 



x^'' exp 



£(5/4), 

r(3/4)' 



fzexp [-2z2(a;-l)2] , 



0, 



CO. 



(73) 
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The crossover of the normahzed distribution P*{x) from the Maxwell-Boltzmann form corresponding to 5 — > — oo to 
the sharp distribution corresponding to (5 = 5 is illustrated in Fig. 




X 

FIG. 12. Plot of the normalized distribution function Ps{x) for 5 = — l)/h^^^ = 
and 5 ( — ). 



(...), -1 (-..-), (— ), 2 (- -), 



D. Critical dynamics 

To study the dynamics in the critical domain define a deviation from the stationary solution by 

V{v,T)^rs{v) [l + y{v,T)]. 

Equation (p7|) becomes 



Now we define a ^[y] by 



Then, from Eq. dTs) we have 



dry (ry, r) = \j (0)^|-P. (,) , Ay r) . 



T[y] = / d7^Vs{v)[yiv,r)Y 
Jo 



^y{v,r) 



(74) 



(75) 



(76) 



(77) 



Thus, J-'[y] has the required properties of a Lyapunov functional for the dynamics of y{ri,T), namely |jTj 
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(78) 



the equality being verified for y = i) only. This implies that for any initial condition the solution to ( [Tq ) evolves in 
time towards 2/(77, t) 0. 

The dynamics for 77(r) is given by 



(79) 



This is not a closed equation so in principle it is necessary first to solve ( |57| ) for the distribution function and then 
calculate rj. However, an estimate can be obtained from (|79| ) using the approximation rf- ^ iv'^s/vDv'^ ■ Then the 
linearized equation for x = 1] ~rj^ obtained from (|79|) for small x is 



drX 



2n*(0) 



x = 0. 



(80) 



Of course, S + 2n^(0)/77, > for (5 > 0. It can be verified from Eqs. ^ and (|2|) that *i(z)/*2(z) > -3 for z < 0, 
s o th at 5 + 2n\{Q)/rj^ > \5\ for S < . This confirms the above stability analysis. Equation ( ^0| ) is consistent with 
(B18) for the ordered phase and (B23) for the normal phase. The finite relaxation time rjg/2nl{0) at the critical point 
is not in conflict with the divergent relaxation time in ( ^l|) since the unit of time is different (i.e., r — \fhs). 



V. DISCUSSION 



In detail application of statistical mechanics methods to the model granular fluid of hard spheres with inelastic 
collisions exposes important differences from normal fluids. First among these is the replacement of the equilibrium 
Gibbs state with the time dependent homogeneous cooling state. In the case of mixtures, the absence of detailed 
balance in collisions leads to a breakdownwn of the usual equipartition theorem for normal fluids. This is interesting 
(e.g., the HCS for a binary mixture has two kinetic temperatures) but is perhaps not too surprising. In the case of a 
single, mechanically different, impurity particle in a one component granular fluid this effect is easily understood as a 
competition between the average impurity-fluid collision rate Vc, responsible for "equilibration," and the cooling rate 
for the fluid ^ constantly changing the reference state. This competition is most severe for conditions such that the 
average impurity-fluid collision rate decreases at constant ^, as occurs when the impurity mass is much larger than 
that of the fluid particle. This requires that the nonlinear dependence of the actual impurity-fluid collision rate on 
impurity mean square velocity is activated to increase the true collision rate. As a consequence, the joint HCS for the 
fluid and impurity is maintained but with a much higher speed for the impurity relative to that for the fluid particles. 
In the limit of infinite impurity mass an extreme breakdown of equipartition occurs with the single impurity particle 
attaining a flnite fraction of the total kinetic energy. 

This peculiar feature distinguishes the conditions of S^/vc — C < 1, where the distribution of energies is similar 
to that for a normal fluid, from ^* > 1 where the distribution is anomalous. A surprising feature of the description 
given here for these two cases is the exact analogy to a second order phase transition in a normal fluid. The order 
parameter is the ratio of impurity to fluid particle mean square velocities, 0^, the conjugate field is h (a measure of 
the mass ratio), and the role of the inverse temperature is the relative cooling rate To summarize the primary 
results obtained here the following comments are offered: 

• The nonlinearity of the impurity-fluid particle collision rate, expressed by the dimensionless friction constant 
7*(0), is essentially the same as that for an impurity in a normal fluid. In the latter case (/) <C 1 for a heavy 
impurity particle and the relevant values are 7*(0) « 7*(0) 1. However, when the background fluid is cooling 
it is necessary that 7*(0) ~ ^* so values of of order 1 are selected when ^* > 1. The details of the mechanism 
by which the host fluid cools is unimportant for this qualitative effect. In fact, even if all collisions are elastic, 
the same two phases would occur if the fluid were cooled by an external thermostat. 

• The thermodynamic analogy originates from an "equation of state" h — h(d>s,£,*), obtained from the "equili- 
bration" condition for the HCS. The phenomenological estimate in Section O and the exact asymptotic kinetic 
theory analysis of Section |^ are essentially the same. The "Gibbs free energy" obtained from integrating the 
equation of state has a Landau-like form near the critical point with critical exponents associated with the 
various first and second derivatives. In particular, the susceptibility diverges, indicating a second order phase 
transition. 



21 



• The approach to the HCS is stable in both phases. The dynamics is governed by a Ginzburg-Landau equation 
defined in terms of the Gibbs free energy. Near the phase transition there is critical slowing, with the character- 
istic relaxation time diverging proportional to the susceptibility. Alternatively this can be viewed as a change 
of time scale from s to t = \/hs. 

• The diffusion coefficient is finite in the normal phase but diverges on approaching the transition. It remains 
divergent in the ordered phase. In terms of the velocity autocorrelation function this is seen to be a divergent 
relaxation time for the decay of correlations, and consequently, the mean square displacement is characterized 
by ballistic rather than diffusive dynamics. 

• The HCS velocity distribution in both phases and in the critical region is obtained from an exact asymptotic 
analysis of the Enskog-Lorentz equation. In the normal phase away from the critical point it is a Maxwellian 
with a temperature different from that of the fluid. In the ordered phase it is a quartic function of the velocity 
centered about a non-zero average speed. The distribution function in the critical region exhibits a continuous 
crossover between these distributions as the cooling rate changes from ^* < 1 to ^* > 1. 

The most direct and controlled observation of the phenomena described here would be via Monte Carlo simulation 
of the Boltzmann-Lorentz equation or molecular dynamics simulation. The qualitative change in the distribution 
function for the ordered phase already has been seen in Monte Carlo simulation. Fig. 6 of Ref. In principle, 
the Monte Carlo simulation could provide access to the longer time behavior associated with critical dynamics and 
diffusion near the critical point. Experimental conditions for real fluids are more difficult to imagine, since a cooling 
medium for the impurity particle is required. However, as noted above, the cooling does not have to be associated 
with inelastic collisions. Thus an impurity particle in a continuously and homogeneously quenched fluid should exhibit 
the same phase transition. 

The extreme breakdown of equipartition discussed in this paper extends to the case of a mixture as well, where a 
mole fraction a;o of impurity particles exists instead of just one impurity particle. In that case a phenomenological 
description similar to that of Sec. |l| shows that the critical value of the control parameter ^* in the limit h —i- at 
finite xo/h is (,* — 1 — {xQ/h){l — ao)/4, so that the HCS of the mixture is always in an ordered state (^* = 0) if 
xo/h > 4/(1 — Oq). The details of this case will be published elsewhere. 
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APPENDIX A: COOLING RATES AND COLLISION FREQUENCY 



The cooling rates ^ and for a- fluid and the impurity particle are deflned by (|^), while the diffusion coefficient in 
Section HI is expressed in terms of a related frequency to d ■ They can be written as 



dt{v\t)) dt{vl{t)) 

so — , 9/, XX ■ 



1 



\t)) 



where the impurity-ffuid particle collision frequency is 

dt' (vo(t)-vo(t + t')) 



(Al) 



(A2) 



The subscript denotes the velocity for the impurity particle and the brackets denote an average over the initial 
ensemble. The time derivatives can be expressed in terms of the generator L for the inelastic hard sphere dynamics 



dtX{t) ^ LX{t), 

N N .... 

i = vo-Vo+^v,.V.-H^r(z,0) + -^^T(z,j). 



N N 



(A3) 
(A4) 
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The binary collision operators for fluid-fluid and fluid-impurity pairs are defined by 

T{i,j) = -a^ j dQ e(-v,, • 5^)(vy • ^)6{r,j - (7){h, - 1), 

T{i, 0) = J dn e(-v,o • a)(v,o • a)5(r,o - o^) (&.o - 1) , 



(A5) 



(A6) 



where cr = a a, cr = aa, and hij and 6io transform the relative velocity for the pairs into their scattered velocities and 
leave the center of mass invariant, 



hj^ij = - (1 + a) (vy • a) ct, hjGij = G 



(A7) 



^ioVio = - (1 + ao) (vio • a-) a-, &ioGio = G^o. 
The various velocities and reduced masses are given by 

v.iO=Vi-vo, Gio = + Movo, /U = ■ , Mo = 



mo 



m + Too 



TO + Too ' 



1 



Vjj=v— Vj, Gij = - (vi + Vj) , 



(A8) 
(A9) 
(AlO) 



= Gio + MoViO, Vo = Gio-MViO, = Gi 



1 



In terms of the binary collision operators the cooling rates and collision frequencies become 



(7v-i)(T(2,iK) + (r(i,oK) 



(All) 



(A12) 



j/o = -N 



(vo-T(I,0)vo) 



(A13) 



In the following it is assumed that terms of relative order can be neglected. Substitution of the definitions for 
T(2,l) and T(1,0) leads to 



T(2, \)vl = -d^ j dn e(-v2i • ?)(v2i • a f&{y2i - a) 



(1 + a) (G21 • ?) - i (1 - a^) (V2i • 5^) 



(A14) 



T(l, 0)^2 = -^24/1 j dfl e(-vio • 5^)(vio • afSino ~ W) {hgw + vo) • ^, 
vo • r(l, 0)vo = -a^2h J dn e(-vio • 5^)(vio • 5^)'<5(rio - o^)vo • a. 



(A15) 



(A16) 



Since these are all two-particle functions the averages in (A_12) and (A13) can be reduced to integrals over the 
two-particle reduced distribution functions /'■^■' and /o^^ defined in terms of the A^-particle distribution function ps as 



(A17) 



f^'^\xo,xi) 



/ dx2 ■ ■ ■ dxNPs (r) . 



(A18) 
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Here V is the volume and Xi denotes a point in the six-dimensional phase space of particle i, i.e Xi <S4> {q;, Vi}. The 
frequencies then become 

^^\na^l-a^)vf J dv^dv; j dnf^^> {wlx^v^i = -ct)Q{^1^ ■ ■ af , (A19) 

x(/ivto+v5).^, (A20) 
1^0 = -2hna'^Vf= J dv^dv^ J dfl /($^^*(v2, v^, rio = -CT)e(vio •ct)(vio • 5^)^Vo • ct. (A21) 



All velocities have been scaled relative to Vf = \/(v^ and 

fi^)=vj^fi^)*, /f (A22) 
The results at this point are still exact. It follows directly from these results that ^ and ujd are manifestly positive. 

1. Neglect of velocity correlations 



If velocity correlations in the reduced distribution functions are neglected on the precoUision hemispheres [£_5 16 
they simplify to 

/(2)*K,v*,r2i = -ct) ^ gf* {vD f* , /f *K,vt,rio = -W) = 50/0 K*)/* K) , (A23) 

where g and 50 ^-i'^ the fluid-fluid and fluid-impurity pair correlation functions for particles at contact. The angular 
integrals can now be performed to give 

f = Irnra^f g (l - a^) J dvldvy* K) f {v^) (A24) 
Co = ^-^^hna^Vfg,^ [{vl'' 1* {v*,^)) h{n*M^))] , (A25) 

y, = ^ (v*) hna'vfgo^ H'l* K')) , (A26) 

where the dimensionless functions 7*(wo^) and ^^(uq^) have been introduced for connection with the discussion in 
Appendix 

7*K^) = J '^vJ/*K)t'oiVo ■ v2i, nlivo^) = J dw*J*{vl)vl^. (A27) 



2. Maximum entropy ensemble 

The HCS distributions are not known exactly, although approximate evaluations suggest they are close to 
Maxwellians. Therefore, to obtain an estimate for the cooling rates and collision frequency the maximum entropy 
(information) ensemble is assumed in this section. This is the Gaussian whose density, momentum, and kinetic energy 
are constrained to have the same values as for the HCS, 

r{v*)=l^^y'e-'^''/^ f*iv*)^(^^y\-'^o'f'~^. (A28) 
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This gives the resuhs 



(A29) 



(A30) 



APPENDIX B: ASYMPTOTIC KINETIC EQUATIONS 

The analysis here is based on the Enskog-Lorentz equation to describe the distribution function for the impurity 
particle. Interest is restricted to the case of small ratio of fluid particle mass to impurity particle mass. To obtain 
an asymptotic form for the kinetic equation, first a Kramers-Moyal expansion is performed to second order in the 
mass ratio. This accounts for the dependence of the coUisional changes on the mass ratio. Subsequently, two different 
expansions are performed for the final asymptotic form depending on the value of a control parameter 

1. Kramers-Moyal Expansion 

The Kramers-Moyal expansion of the Enskog-Lorentz equation has been obtained in Appendix A of reference [|j . 
The result is 

d 19^ 

9i/o(vo,i) = ^ [hvo'j{vQ)fo{vo,t)] + -- — - — {/i^ [ni{vo)S^j 

dvo 2 dvoidvoj 



+n2{vo) [ voiVoj - ^Sijvl 



The friction 7(wo) and noise ni{vo) , n2{vo) are 



3Vr 



3h 



5h 



where Vc = ^hmrcr'^go (v) i s the characteristic impurity collision frequency introduced in 
have been defined in (A27), and 
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16 (v* 



I * ^ \2 1 *2 



The dimensionless variables are 



(Bl) 
(B2) 

Also 7*(0) and nl{(j)) 
(B3) 

(B4) 



The analysis of the deterministic limit in the text makes use of the property ^*{4>) > 7*(0) = 1. To prove this, first 
perform the angle integrations to get 



7*(<^) = 1 + 



> 1 + 



1 



5{v*) 
1 

5{v*) 



(w*-^) -n<f) dv*f*{v*)v* 4 -I- 



V0 



(w* — 47r0 / dv*f*[v*)v* 



(B5) 



where the inequality results from 4 > (4 + (1 — x)^ for a; < 1. Next, writing out the contribution from (w* ^) 
explicitly gives the desired result 



7*(0)>1 + 



5{v*) 



_dv*f*iv*)v* > 7*(0) = 1. 



(B6) 
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Analogously, it is possible to prove that 7*'(0) > 

1 



7*'(0) 



> 



5{v*) 
1 

47r 



*2 



5 (v*) 



(w*-^) -47r / dv*f*(v*)v* 
Jo 

dv*f*{v*)v* > 0. 



V0 



(B7) 



The remaining analysis of the text and below does not require the explicit forms for 7* (</>), (</>), and 71-2 Ho wever , 
for the illustrations in the graphs an excellent approximation is obtained using the maximum entropy ensemble ( A28| ) 
for the fluid; no assumption is required regarding the impurity particle distribution. The resulting integrals can be 
performed with the results 



7*W = ^(l + 3</')e 



-30/2 



1 /27r , 



902) erf v/3072 , 



(B8) 



/2F 



902 erf K/30/2 , 



(B9) 



5 5 /27r / \ 

"2(<^) = 480^ + 20 + 302) + (1 - 30 + 902 + 90^) erf (^3072) • (BIO) 

In the following only solutions that depend on the magnitude of vq are considered. Since the order parameter is 
the average of 0, 

,,*2 



0(t) = J d^rov*o'fo{vo,t), (Bll) 

it is appropriate to change variables from vo to 0. In addition, the dimensionless time scale of (|^) is introduced. This 
is accomplished by defining the new distribution function P(0, s) by 

P(0,s) = A7Tfoivo,t)vl^ = 2^«3^i/2/o(z;o,t), (B12) 



The Kramers-Moyal equation ( |Bl| ) becomes for P(0, s) 

a,P(0, = ^ {-^ [-r + T*!-^)] - (1 - l^^^ hnli^) + ^^^^hn*icf>)^ P(0, s) 

+0{h^). (B14) 

The deterministic limit, ft. = 0, is described in the text. In the following an outline of the fluctuations about this 
determistic limit is given. 



2. Expansion around 0qj, 

The effects of finite h represent "noise" which broadens the width of the initial delta function as the system evolves. 
To include such effects consider solutions of the form 
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P(0,s,/i) = h-PVi- 



hP 



(B15) 



such that the hmit hm/i^o 'P{rj^s^ h) — V{rj,s) is finite and independent of h. The choice of reference state given by 
( pT| ) imphcs initial condi tions that do not deviate too much from the stationary state. To find such solutions, define 
a change of variables in (B14) by 



= 4s + h^V, Pi^: s.. h) = h-PV (ry, s, h) . 
In the ordered phase a non-trivial equation for V is obtained with the choice p — 1/2, 



5,7'(r,,s,0) 



d r - 



dr] \ 



?/0Os7*'(0Os) + (t^o 



2 — 4— — ' 



-I 



P(ry,5,0), r>l- 



(B16) 



(B17) 



where it has been recognized that i^* — l*{<t>s) = in this phase. The average value of -q obeys the equation 

5,77 (,s) = -0037*' J77(s), (B18) 
for ^* > 1. Stability is assured by 7*'(0os) ^ 0, which 



which is the linearized form of the deterministic dynamics (BO 



limit. The stationary solution for the distribution function is obtained from 



is seen to be the case using (B7). The stationary solution to (B18) is 77^ = so there are no corrections to in this 

(B19) 



2 — 4_ _ 



d 



= -0os7*'(</'os)^n, 



whose solution is 
or, in terms of 



1 



B I 



7*' (00s) 



2 — 



4- 



2B(0oJ/i^ 



(B20) 



(B21) 



In the normal phase — and the change of variables in ( B16 ) becomes (f> = hPrj. A non-trivial equation for V 
is obtained with p — I, 



The average value of 77 now obeys the equation 
The stationary solution to this equation is 



Vs 



which gives the leading finite contribution to 0s as /i ^ 0. The stationary distribution function is 



3 f 3ri 



and the corresponding distribution in terms of is 



^^s(0) 



30 



1/2 



1/2 



,-30/2h>7. 



(B22) 
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